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Abstract 

An effective Hamiltonian is derived for two coupled three- Josephson-j unction (3JJ) qubits. This 
is not quite trivial, for the customary "free" 3JJ Hamiltonian is written in the limit of zero induc- 
tance L. Neglecting the self- flux is already dubious for one qubit when it comes to readout, and 
becomes untenable when discussing inductive coupling. First, inductance effects are analyzed for 
a single qubit. For small L, the self-flux is a "fast variable" which can be eliminated adiabatically. 
However, the commonly used junction phases are not appropriate "slow variables", and instead one 
introduces degrees of freedom which are decoupled from the loop current to leading order. In the 
quantum case, the zero-point fluctuations (LC oscillations) in the loop current diverge as L — > 0. 
While their effect thus formally dominates over the classical self-flux, it merely renormalizes the 
Josephson couplings of the effective (two-phase) theory. 

In the coupled case, the strong zero-point fluctuations render the full (six-phase) wave function 
significantly entangled in leading order. However, in going to the four-phase theory, this un- 
controllable entanglement is integrated out completely, leaving a computationally usable mutual- 
inductance term of the expected form as the effective interaction. 

PACS numbers: 85.25.Cp, 85.25.Dq 
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I. INTRODUCTION 



A commonly considered flux qubit consists of a superconducting loop with three Joseph- 
son junctions — a 3JJ qubit One readily writes down its Hamiltonian, 1 

3 
i=l 

in units with H = 1, and with the external flux bias $ x given in phase units, 

X = 2e$ x . (1.2) 

As long as one neglects gate capacitors etc., the Ci are simply the junction capacitances: for 
finite inductance L, the redistribution of charges following a tunneling event, often accounted 
for by an effective capacitance matrix, is not instantaneous, i.e., all three junction charges 
and phases are independent dynamical variables. 2 

The essence of the 3JJ design is to introduce bistability without relying on magnetic 
energy, so that the SQUID loop can be kept small. For L — > 0, the last term in (jl.lj) mainly 
implements the constraint 03 = X — <f>i — 02- This leads to an effective two-phase theory, 

#3JJ = ^Q^C^Q - Ei COS 01 - £ 2 COS0 2 - £3COs(0 x -0i-0 2 ) , (1.3) 

fC x + C 3 C 3 \ , > 

\ c 3 c 2 + c 3 J ' li - 4j 

with Q = (Qi, Q 2 ) T (cf. (4) in bearing in mind that 0i and 2 have the opposite relative 
sign there). 

When scaling to circuits with more than one qubit, the question arises how to generalize 
the effective Hamiltonian (jl.3|) . Unfortunately, repeating the above reduction merely yields 
uncoupled 3JJ Hamiltonians. Clearly, by neglecting self-inductances entirely one has missed 
the effect of mutual inductance (which can never be larger), responsible for qubit-qubit 
interaction and ultimately entanglement. This can be turned into an advantage: most of 
the analysis can be done by studying one qubit to higher order in L than usual, after which 
two qubits require only a straightforward generalization. 

This is taken up in Section |TI] on the level of classical Kirchhoff (circuit) equations. 
Section llTTl is devoted to the quantum case. Rather than merely confirming the classical result 
or resolving operator-ordering ambiguities, one finds an additional physical effect: zero-point 
fluctuations in the loop current wash out (renormalize) the effective Josephson couplings. 
The results are verified numerically, and critically compared to a previous attempt |2j. 
Section II VI contains the generalization to two qubits; it is seen that the four-phase wave 
function is entangled only in considerably higher order in the inductance than the six-phase 
one. This integrating-out of uncontrollable entanglement may indicate that the effective 



1 As always, <fii in Ijl.lD are properly gauge-invariant "phases" , defined as line integrals along the SQUID 
loop. Thus, there is no reason for H to be 27r-periodic even though the Josephson term is, and indeed 
(pi's differing by multiples of 2tt correspond to physically distinguishable magnetic field configurations. 

2 A nontrivial capacitance matrix will be essential in Note also that, for this pure series circuit, 
H does not depend on how the inductance is distributed along the loop. 
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theory is not merely computationally convenient, but also physically appropriate. Some 
concluding remarks are made in Section IVl 

For small L, self-fluxes vanish but persistent currents remain finite, so that the latter 
are often more convenient. When studying dynamics etc., one needs the current operator 
(as opposed to, say, the ground-state expectation ||), which is not readily available in 
the growing literature on the 3JJ qubit. Preliminary to the L-expansion proper, it will 
now be derived quantum-mechanically in both the three- and two-phase theories, which has 
some independent interest. In fact, one merely needs to add the capacitive and Josephson 
contributions. Starting from 1)1.1)1 . one has 3 



h = -hi sin ^ - Qi 



with I r j = 2eEj and 



Then, 



Qi 



2e d 



Qi = [H,iQi 



2eL 



- Id sin , 



2eL 



(1.5) 
(1.6) 

(1.7) 
(1.8) 



That is, h = I independent of i as expected; contrast, e.g., arbitrarily picking the Josephson 
current through one of the junctions. 

For L — > 0, p. 8)1 tends to jj and a separate derivation from ()1.3)1 is needed. In this case, 
however, the canonical variable Q± is not simply the charge on capacitor 1 (indeed, two Q's 
have to account for three C's); to find the latter, we set C _1 = (rq) and evaluate 



ii = i[±pQ?+rQiQ 2 , 0i] = 2e(pQ 1 + rQ s 



(1.9) 



Hence, 



/{ 0) = 
1 2e y 



J cl sin 0i 

(Cip-l)/ci sin 0i + C x rl c2 sin0 2 - Ci(p+r)/ c3 sin(0 x -0i-0 2 ) 
-C 



hi ■ , Ic2 ■ , . Ic3 . i 

— sin 0i + — sm 02 + — sm( 

Gi O2 O3 



1.10) 



an appealing and manifestly symmetric form I\ 



(o) 



r(0) 



involving the loop's series capaci- 



tance C 1 = C x + C 2 + C 3 . Interestingly, the junction area cancels from p. 10)1 : if, 
as is often assumed, I ci oc Cj oc Kj, then the three sines have equal prefactors. 



3 When choosing conventions, one needs a minus sign somewhere for the final 1)1.8(1 with ((1.2JI to come out 
as $ - <f> x = LI. 
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II. CLASSICAL ANALYSIS 



The Hamilton equations 



2e 



OH 



Qi = -2e 



dH 

d<pi 



(2.1) 



yield the classical dynamics of the system (jl.lj) . 4 The last term in H represents a steep and 
narrow well, so one expects 

= 01 + 02 + 03 ~ 0x (2.2) 

to be small and rapidly oscillating, while the other variables can have excursions of order one 
but move comparatively slowly — ideal for adiabatic elimination of 0. However, from ()2.1|) 
one finds, e.g., Ci0i = —(j)/L — ZeldSmtpi, so that determination of the 0i-dynamics to 
O(L) apparently involves the 0-dynamics to 0(L 2 ), rendering the L-expansion quite tedious. 

To see what went wrong, consider the rough analogy of a two-dimensional electron gas, 
where excursions away from the x-y plane carry a large penalty in energy. It is intuitively 
clear that the fast oscillations in the corresponding potential well occur predominantly in 
the z direction, orthogonal to the easy plane. While the latter plane is uniquely defined 
by the potential, in our case "orthogonal" involves the anisotropic capacitance matrix as 
well (which thus can be said to induce a metric). It is not difficult to scale the phases such 
that the charging term is oc Q T Q (in terms of a three-vector Q; in the quantum case, this 
yields a charging term cx —V 2 ), upon which the proper coordinates are found by rotation. 
However, full orthonormalization turns out to be overkill, and the resulting tedious formulas 
are hard to interpret; crucial is only that the slow coordinates x> @ are constant along the 
fast direction 



n, <f>2, 



(C 1 ,C 2 X ,C 3 x ). This is readily achieved; the easiest seems 5 
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The Kirchhoff equations become 

CiX 
2e 

c 2 e 

2e 



cr 1 cvV 



J c3 sin 03 - J c i sin 0i , 
7 c3 sin 03 - J c2 sin 2 , 




(2.4) 



(2.5) 
(2.6) 



4 Since superconductivity is not a "classical" phenomenon, the term is a bit tenuous. Indeed, if H is written 
in terms of phases (fluxes), the inductive (Josephson) terms would involve h in SI units. Of course, the 
dynamics (|2.1|l do emerge from the full quantum^ theory in the limit of large capacitances; cf. Section Ivl 

5 The coordinates 0i, 82 (plus I m oc 4>) in (13) of [2j satisfy the same criterion. They correspond to choosing 
two zero entries in (12.4(1 rather than (|2.3(l . and are physically equivalent to x, ®- 
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where the 0,'s on the rhs are simply shorthand for the linear combinations in (12. 4|) . Indeed, 
the terms ~ L~ l have cancelled in ()2.5|) and ()2.6|) . For an interpretation, note that x (X 
C\Qi — C3Q3 describes the charging of the island between junctions 1 and 3; unlike the 
charging of, e.g., capacitor 1, \ thus has no net direct contribution from the loop current, 
and is affected by it only indirectly [in O(L)] through the flux-quantization condition. 

In these variables, the L-expansion to the required order is almost trivial. In O(L ), one 
simply sets = in (J2.5j) and (J2.6|) . which is readily verified to yield dynamics equivalent 
to (jl.3j) . In O(L), can be neglected on the rhs of (J2.7|) since it does not occur with a large 
coefficient. The leading adiabatic solution is then 



2eLlV>(x,e) 



(2.8) 



or explicitly l'°>(Xi 0, (f>=0) in (jl.lOj) . Subsequently, ()2.8|) can be substituted into the rhs's 
of (|2.5|) and (|2.6|) . yielding a self-contained dynamical system with two degrees of freedom. 

In preparation for the quantum analysis, it is instructive to also consider a canonical 
formulation. The Lagrangian reads 



1=1 

c 

8^ 



8e 2 



V 



(2.9) 
(2.10) 



in which one is to substitute ()2.8|) for 0. The transformation (|2.3|) has achieved that ()2.10|) 
does not contain x<fi or 6*0. Such terms, which do occur in C((f>i, 02, 0), would in 0(L) result 
in an extremely unpleasant phase- dependent capacitance matrix in the effective Lagrangian 
(or Hamiltonian) . However, in ()2.10|) with ()2.8|) . 2 (x, #) is 0(L 2 ), hence negligible for the 
(x, ^-dynamics. The effective Hamiltonian follows from Q x = 2edC/dx etc., as 



H(x,9;Q x ,Q e ) 



Q X X + QeO 



2e 

1 

C~i 



9i 
'c x 



2C. 



1 

C~2 



9i 
c 2 



The potential V merits closer inspection: to O(L) 



v(x,e) = u J (x,9,cj ) =o) + 



4>=o 



8e 2 L 



+ V 



(2.11) 
(2.12) 



(2.13) 



Using ()2.4|) to perform the differentiation in the second term explicitly at constant x, d, one 
finds fyC/j| = -J(°)(x,^)/2e. 6 Substituting (22J) for </>(x,6), this means that this term 



By the same token, one has Z(°) = -2edH /dcpx\ Xie for H as in H2.15J1 . This seemingly obvious relation 
is simply not true for dH^j / d4>^\ t j J1 ^ 2 in even though I = —2edH/d(j) K \^ does hold between Ijl.lfl 

and l|1.8ll in the three-phase theory. For coupled qubits, one likewise has dH/dB x = —A a I a — A^Ii, in 
l|4.4ll . (|4.5|l below. 




(a) (b) 

FIG. 1: The zeroth-order potential Vq (a) and the first-order one V (b), for a = 0.8 and <^ x = n. 
In (b), we also use e 2 LE\ = 0.3. The contours correspond to Vt^/Ei = —1.42, —1.3, —1.2 
(= a— 2), —1, —0.8 (= —a), —0.5, 0.5, 1.5, and 2.5, the latter encircling the maximum at 
X = ^ 2 = tt. The well V (±<p* ,±f f )/E 1 = -a - l/2a = -1.425 at £* = arccos(l/2a) ps 0.8957 
is shifted to V(±g f ,±g')/Ei « -1.7983 at ps 0.9602. While the shift in well location thus is 
appreciable, just evaluating the persistent current in the unshifted one, i.e. V(±4>* ,±(j)*)/Ei = 
—a — l/2a — e 2 LEi(2—l/2a 2 ) —1.7906, gives a good estimate of the well depth. 



inverts the sign of the third one, yielding 
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(2.15) 

'Josephson" 
2,0), where 



Here, Hq plays the role of H^jj in (%, 9) coordinates. The combination of 
and "inductive" corrections in ()2.14|) would not have occurred in, say, (0i, 
thus an awkward potential term would have compensated an awkward charging term [cf. 
below (j2.1(Jj) ]. The sign-flipping of the magnetic term has also been observed around (2) 
in Ref. |4j, though not in terms of the detailed current operator l' ' in (jl.lUJI . Thus, the 
semiclassical analysis in p| seems a valid approximation for stationary currents in the low- 
lying qubit states. 

Since H is invariant, one can here and in the following increase the similarity to the 
conventional treatment by introducing 4>i = C^C^+C^x — 9jC\ + 4> X /Ci} and </> 2 = 

C{— x/C 2 + (C± 1 +Cz 1 )9 + X /C 2 }, so that formally Hq{4>) = if 3JJ (0). Since comparison 
with ()2.4|) and ()2.8j) shows that 0i i2 do not coincide with the junction phases to O(L), this 

formulation is not very useful in derivations. However, in terms of 0, the double periodicity 
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of V and any additional symmetries are expressed conveniently, which may be effective 
numerically [cf. below ([3.28)1 ] . Figure ^ shows both the zeroth-order potential Vq and the 
corrected V of ()2.13|) . for standard parameters E\ = E 2 , C\ = C2, and a = E^/E\ = 
C3/C1 = 0.8. In this case Vo(0i,02) = Vo(02,0i), which is preserved in V; since we have 
chosen a degenerate bias X = ir, one also has V(o)(0i, ^2) = V(o)(—<Pi,—<pi)- The small 
self-flux has a clear physical effect: (I^) 2 = V^(/ (0) ) 2 = for T e {(0, 0), (0, tt), (tt, 0)} by 

symmetry, so both the location and the height of all saddle-points are unchanged. 7 Since 
clearly — ^L(I^) 2 < in the wells, both the "easy" (intra-) and "hard" (inter-cell) barriers 
are increased, and this term's suppression of the tunneling amplitude may not always be 
negligible. Finally, the lines 0i + 02 — 7r (mod 2ir), which are straight equipotentials 8 
through the inter-cell saddle-points of Vq, cease to play this role in V. 



III. QUANTUM ANALYSIS 

A. Expansion of the Schrodinger equation 

In quantum mechanics, one can advantageously use the same variables ()2.3j) . avoiding 
<9 X <9^, and ded^ in ()3.2j) below. The large inductive term in (jl.lj) for H now causes the system 
to remain in its ground state with respect to 0. To leading order this is a harmonic-oscillator 
ground state ip ~ e~ A * , and one readily finds 

This unfortunately means that formally <p ~ L l l 4 (d^ ~ L" 1 / 4 ) and strikingly / ~ L~ 3 / 4 . 
Hence, the expansion to O(L) will turn out to be sixth-order while it was first-order in 
Section |TIJ and keeping the calculation organized is essential. Since it suffices to study the 
time-independent problem, the expansion will be analogous to, e.g., fast- variable elimination 
in (Fokker-Planck) diffusion operators p|. 

Using (jl.fij) and ()2.3j) in (jl.lj) . in the eigenvalue equation (H — E)ip = we expand 

H = - 2 ^ d l + ^Z + H * + H ^ + H ^ 2 + --- > (3-2) 
E = E zp + E + Ex + ■ ■ ■ , (3.3) 
i> = ^ + ^2 + ■ ■ ■ , (3.4) 

with the zero-point energy E zp ~ L~ x l 2 and H as in (|2.15|) . with Ei,ipi ~ L l//4 , and where 
we have been able to omit an L -1 / 4 (L 1 ^ 4 ) term in E (ip) from the outset because H does 
not have an L~ 1//4 term. One has H n = (l/n\)d^\ Q Uj for n > 1, and since the calculation 
below (|2.13|) is just as valid in quantum as in classical mechanics, Hi = — /^(x, 9)/2e. 



7 The same holds for the maximum V^(tt, tt)/ ' E\ = 2+a, but this should have few physical consequences. 

8 This may be a new observation, and holds for general tfi x even though the saddle-points themselves will 
move along this line. 
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In O^L- 1 / 2 ), one has 



2p 2 

Q °<t> + 



8e 2 L Ezp 



^o = 



zp 



(3.5) 
(3.6) 



the obvious answer for zero-point fluctuations of the loop current. 

One can use ()3.6|) to factorize the operator in brackets in ()3.5|) . In O(L ), this leads to 



2^ 
' C 



(3.7) 



where one can commute e Xf ^ 2 through Hq—Eq. Operating with Jd0e x ^ 2 , the first term 
vanishes, so that one obtains the solvability condition 



(H -Wo = 



(3.8) 



equivalent to the standard 3JJ theory (jl.3|) . Substitution back into (J3.7|) renders the latter 
identical to the leading order, so that 



^ 2 (x,M)=V4(x,#)e- A02 • 
Proceeding to (^(L 1 / 4 ), one has 



1(0) 
'~2eT 



(3.9) 



(3.10) 



for a solvability condition 



E 1 = (3.11) 
because the term oc I^cfi cancels by parity. This term does however contribute to ip3 itself, 9 



^s(x,M) 



16e 3 A 



(3.12) 



where the first term in square brackets is readily verified to shift the center of the Gaussian 
in (|3.6|) to 10 (<fi) x o i n accordance with (|2.8|) . 

In 0{L x l 2 \ one obtains the first nontrivial solvability condition, 11 



(# -£ )^+(|l-£ 2 )^ = 



(3.13) 



9 Only solutions bounded in <fi are acceptable. 

10 Properly, (<^) Xo g = (<M(x~Xo)<5($ — $o))/(£(x~Xo)c>(#— #0)) is a conditional quantum-mechanical expec- 
tation. The vector version below l|4.10|l is analogous. 

11 Of course, one can operate on (|3.13|) with f d%d# (ip' )* , yielding E 2 = (iP'q\(H 2 / 4:X)\ip' } / (ip' Q \^' ) , upon 
which ip' 2 is determined up to a multiple of tp' Q (a trivial change in normalization). However, our focus is 
on deriving the effective theory, not solving it. 
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where the prefactor is obtained as f d(p (p 2 e~ 2X ^ 2 = (4A) _1 /d0e~ 2A ^ 2 . While we are only 
expanding the linear Schrodinger equation, it apparently "knows" that the proper statistical 
weight is \ipo\ 2 not ipQ. The correction to ip follows as 



CH 2 ip' 2 , , m 



16e 2 A 

In (9(L 3//4 ), the solvability condition reads 

(H -E )i/>' 3 -E 3 iP' = 



-\<j> 2 



(3.14) 



(3.15) 



apparently posing the obstacle that both ip' 3 and E 3 are unknown in this order. However, 
one can operate with J d\d9 {ip' Q )* and use ()3.8j) . yielding 



and "03 oc ip' , so that 



En = 



(3.16) 



(3.17) 



is a convenient choice of normalization. Determination of ^5 presents no problems, but is 
not required. 

Finally, one arrives at O(L), finding a solvability condition 

(H -E )^ - e 2 L(H -E )H 2 ^' - ^(/ {0) )Vo + (f^V + (j^H 4 -E^' = . 

J J (3-18) 

What remains is to interpret this in terms of an effective two-phase Schrodinger equation. 



B. Effective theory 

While odd powers in L 1 / 4 were important in the above, most clearly through J(°) in (|3.12|) 
implementing the quantum counterpart of tj2.8j) . they cancel in final results such as (j3.11|) . 
(I3.16|) by </>-parity. What emerges is a two-phase theory involving an expansion in terms of 
in which the zero-point motion ()3.5|) is absent, the leading order being simply 

(H -E )^ = 0, (3.19) 
H = H , E = E , ^o = ^o- (3.20) 

Subsequently, (j3.13J) is recognized as the first-order part of 

[H 0+ H 1+ (£ +£i + • • • )](^o+^i+ • • • ) = , (3.21) 

if one sets 

= E 1 = E 2 , ^1 = ^2- (3-22) 

For the next and final relevant order, care is_ needed in identifying the effective wave 
function. Since probability should be conserved, ip must obey \i/i(x, &)\ 2 = f d</> \if)(x, @, 0)| 2 , 
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in which ^3 cancels to 0{L) by (gH and (l3~T7f) . Combining (jS3j), (|3~K|) . (Q, and (f3~Tl) . 

one has 



|^| 2 + 2 Re[(^)> 2 ] + |^ 2 | 2 + 2 Re[($>) V 4 ] " 2e 2 L# 2 K| 2 + 0(L 3 / 2 ) = $ +$ l+ $ 2+ . . . | 2 , 

(3.23) 

which is satisfied by 

^ 2 = - e 2 L# 2 ^o . (3.24) 
Substituting this into (J3.18j) . the latter is identified as the second order 12 of (J3.21|) . with 13 

E 2 = E,, H 2 = l ^H,-^f . (3.25) 

Combining the above, one sees that the effective quantum Hamiltonian H again has the 
form (j2.14|) . but with renormalized Josephson couplings 



^^-fVp + fg!). ( , 26) 

To this order, the latter is exactly what one would expect from a Gaussian averaging 
U 3 (X,9) = Jd(f)e- 2X * U 3 ( X} e } (j))/ Jd(f)e- 2 ^ : the zero-point fluctuations in behave as 
uncertainty in X in (|2.15|) . which washes out the Josephson potential. The systematic ex- 
pansion is not prejudiced about "Josephson" vs "inductive" contributions to the energy, and 
instead directly provides a term — ^L(I^) 2 to H with the correct sign. 



C. Numerics 



To test the above analysis numerically, introduce the dimensionless inductance (3 = 
27rLJ c /$o = 4e 2 LEj ($ is the flux quantum) and g = Ej/Eq (with the charging en- 
ergy Ec = e 2 /2C\), and measure all energies in units Ej. Evaluation of the L-expansion 
is only meaningful for (3 < 1, and j3 should be varied by at least an order of magnitude. 
These imply that for some of our parameter values, the energy levels of the full and effective 
theories will be very close. Hence, we focus on the simplest junction systems (to which the 
above for the 3JJ qubit is readily adapted), stressing accuracy rather than generality. 

First, consider the rf-SQUID, H = — (4/g)<9| — cos(0+0 x ) + 2 /2/5, for which the effec- 
tive H = a/2//% — (l — ^/f3/2g+(3/Ag) cos0 x — ^(3sm 2 (f) x is zero-dimensional, i.e., a closed 
expression. 14 To evaluate H, we use that <fi has an effective range ~ A -1 / 2 [cf. (|3.1j) ]. by 
truncating the potential with hard walls at <fi = ±m^/(3/g. Discretizing the resulting Dirich- 
let boundary value problem, the ground-state energy E is found by direct diagonalization, 
checking for convergence both with respect to the grid spacing A and m = 6, 7, 8,.... 15 



12 The occurrence of H\ipi in this order represents the effect of H\ in second-order perturbation theory, not 
a new term in H . 

13 Again, one can solve E2 = [(V'ol-f^lV'o) + (ipi\Eo— -ffo|'^i)]/(?/'o|V'o)) cf. note^2 However, in practice direct 
numerical solution of H — though no more accurate — seems preferable over a perturbative approach. 

1 For convenient comparison, E zp (usually an irrelevant constant) is included in H. 

15 Increasing m too fast, e.g. successive doubling, wastes grid points on the wave function's large-|0| tail. 
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40 


7T 


0.32 


0.471441 


559017 


463997 


40 


7T 


0.08 


1.097919 


1.118034 


1.097038 


40 


7T 


0.02 


2.231116 


2.236068 


2.231006 


40 


7T/4 


0.32 


— 1 057313 


— 1 076987 


—1 059445 


40 


7T/4 


0.08 


-0 503724 


-0 517970 


-0 503857 


40 


vr/4 


0.02 


0.608463 


0.600064 


0.608456 


1 


vr/4 


0.32 


2.998577 


2.870307 


2.996471 


1 


vr/4 


0.08 


6.476253 


6.405841 


6.476015 


1 


vr/4 


0.02 


13.513831 


13.476908 


13.513799 



TABLE I: The ground-state energy E of a symmetric dc-SQUID, compared to the standard esti- 
mate Eq which ignores inductance effects entirely, and to the prediction of the effective theory E, 
accounting for renormalization of the Josephson coupling and for self-flux. Special attention has 
been paid at g = 40, X = 7r/4, /3 = 0.02 to obtain the very small E — E reliably. For X = tt, the 
zeroth-order problem reduces to a free rotor, hence Eq = 2/^/]3g. 



Using the known error ~ A 2 of lowest-order discretization of <9|, one can perform a (highly 
effective) Richardson extrapolation .E(AjO) ~ [4i?(A) — E(2A)\/3. The resulting data (not 
shown) confirm that the effective theory has an error 0((3 3 ^ 2 ). 

For a system with a non-trivial effective theory, consider the symmetric dc-SQUID, H = 
— (4/#)(c^i+c^ 2 ) — cos 0i — cos 02 + (0i+0 2 — x ) 2 /2/3. With fast and slow variables = 
01+02— X , X = — 02) respectively, this becomes 

2 

g "*-- x 2 J™"" 1 2/3 



if=--9j--a£-2cos( )cosx+^ , (3.27) 

$7 <f x "V" 4V2'32s;-~V2y"~ /v 2""V2 




2 n2 1 P , "\ /? . 0/ "x \ 2 



if = -^--d 2 -2 1-tW-+^- cos cosx-^-sin 2 cos 2 x • (3.28) 



In both (|H.27jl and (|H.28j) . 27r-periodic boundary conditions are imposed on \. This is not 
the only possible choice, and corresponds to a zero offset charge in the two arms of the 
SQUID;_for large g, the effect of this boundary condition will be small. Observing that both 
H and H are even in x, the even sector (containing the ground state 16 ) can be treated as a 
Neumann boundary value problem on < x < 7r - Further details are as for the rf-SQUID. 17 
Some results are shown in Table HI together with the standard zeroth-order prediction Eq 
which results by dropping all correction terms in f)3.28j) . In all cases, reducing (3 by a 
factor 4 reduces the error E — E by at least a factor ~8, as predicted. We note without 
explanation that for g = 40, X = 7r/4, this factor seems to be rather ~16, as if the next 



The ground state tends to be found with the highest accuracy for a given A, but other states are readily 
found and compared as well. 

The used matrix representation [6j makes it convenient to allocate only half of a symmetric matrix in a 
way that is transparent to the rest of the program, as long as the latter never accesses matrix elements 
above the main diagonal. 
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0(j3 3 / 2 ) correction were very small numerically or cancelled altogether. In contrast, E — E 
only goes down by a factor ~2, except for <f) x = it, in which case the effective theory has 
a zero Josephson coupling so that its 0(\/j3) renormalization does not matter. Still, even 
in the latter case the self-flux term is nontrivial, so that an 0{f3) error renders Eq less 
accurate than E. The last three rows emphasize that, in spite of this paper's title, our 
//-expansion is not at all limited to the "flux-qubit" regime g ^> 1. This is useful, since 
small-L, small-C dc-SQUIDs are commonly used as tunable compoundjunctions, and 3JJ 
qubits in an intermediate-^ regime have also been considered recently [7(. 

In summary, the transparent effective Hamiltonian H can be used with excellent accuracy 
also in many cases where the zeroth-order H would be unsatisfactory. In the actual com- 
putations, the advantages in speed and storage requirements are evident; these will persist 
even when less naive algorithms are employed in the study of more complex devices. 



D. Comparison 

Inductance effects in flux qubits have previously been studied in Ref. jH, 18 with the 
examples of loops containing one through three junctions. For one junction, the rf-SQUID, 
the two treatments are still largely equivalent, and agree on the 0(\/L) contribution to 
the energy [due to the second term in ()3.26|) ]; this presumably is the main reason for the 
improved agreement with numerics in their Fig. 1(b). In O(L), the last term in (5) gives 
the self-flux contribution — |L(/^) 2 ; however, in the same order the renormalization effect 
[third term in ()3.26|) ] is missing, because in the transition from (3) to (4), the Josephson 
potential was expanded to second instead of to the required fourth order in /. 

For loops with multiple junctions, several problems arise in the treatment of 0. First, 
the physical role of the loop's series capacitance is not identified, 19 as is most clearly seen 
in the third term of (9), written in terms of the parallel capacitance 2C (the appropriate 
effective mass for the G coordinate) while this term can be naturally expressed as |(C/2)L 2 / 2 
[cf. (J3.27J) ]. Observing that the series capacitance of the 3JJ qubit is [a/(l+2a)]C would 
have uncovered the missing factor | in the third term of (14), which presently carries through 
to underestimating u below (16) by a factor v2 — which is not detected in Fig. 3(b) because 
the latter does not feature three-phase numerics. 

Second, the renormalization effect is not estimated correctly. The last term of (10) is 
neglected because it involves L 2 , which is small, without checking whether I 2 could possibly 
be large. In fact, this term's counterpart in (4) yields the dominant correction in (5). For 
three junctions, the corresponding term is simply omitted from (15). Note that Fig. 2(b) 
is given on such a scale, and for such extremely small (3, that the apparent agreement with 
numerics only verifies the zeroth-order result which never was in doubt. 

Third and most serious, the analysis is conceptually incorrect in averaging the dynamics of 



In that paper's opening paragraph, 27rJ c /$o defines LJ 1 not Lj. The order estimates below (1) give the 
number of grid points, while the "computational time" grows faster. Equation (3) and the one right above 
have a sign problem when compared to (2) and Fig. 1(a). A constant 2Ej is lost between (6)-(7) and (9). 
In the definition of m below (14), the factor 2 in the second term is spurious. In the expression for hujQ 
below (16), E c /Ej should read E c Ej. The figures are irreproducible since Ej/E c is not given. 
In general, not taking independent capacitances and Josephson couplings makes it more difficult to con- 
ceptually separate the effects of the two. 
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the fast variable over the dynamics of the slow ones [first below (9)] instead of the other way 
round. In the final results (11) and (16), this leads (in our notation) to a perturbative energy 
correction — ^L(I^) 2 instead of the proper — |L((/^) 2 ) (cf. note IT3*j) . This is inconsistent 
with the classical limit of Section|n]and underestimates the magnitude of the self-flux effect- 
most dramatically for the two lowest eigenstates of the degenerately biased 3JJ qubit (cf. our 
Fig. [I}, where it would incorrectly predict this effect to vanish entirely. See also the remark 
on the coupled case in Section M below. 



IV. TWO COUPLED QUBITS 



Consider two SQUIDs a and b side by side. The directions of the fluxes have to be chosen 
consistently in each loop w.r.t. the external field, and their magnitude is given as 

$-$ x = LJ, (4.1) 

where <fr = ($ a , ${,) T etc., and where 

so that the antiferromagnetic coupling is characterized by a positive M. We can assume a 
homogeneous external field, so that the applied fluxes <& x = (A a , A b ) T B x are proportional to 
the loop areas A a>b . The magnetic energy can now be written as 20 

H M = |i" T L/ (4.3) 

($ a -A a B x ) 2 (<& b -A b B x ) 2 M($ a -A a B x )($ b -A b B x ) 
2L a (l-P) + 2L b (l-k 2 ) + L a L b {l-k 2 ) ' 1 ' } 

where k 2 = M 2 /L a L b <C 1 will never (need to) be assumed. The factors (1— A; 2 ) -1 S\ are 
often ignored, but should not be surprising given the nontrivial effective capacitances one is 
used to seeing in charge-qubit Hamiltonians. Physically, they arise since, e.g., $ a — A a B x is 
not simply the self-flux of loop a, but includes a mutual contribution. 21 



Expanding the bilinear form in 14.3(1 . the aniiferromagnetic interaction is seen to involve a term ~MI a Ib, 
in contrast to the positive last term in (|4.4(l . This is counterintuitive, but a simple consequence of 
2x2 matrix inversion. It may be clarifying to compare (anti)ferromagnetic configurations at constant 
absolute fluxes in the loops. Due to the mutual inductance, antiparallel fluxes can be generated with 
smaller |/|'s, reducing the value of ^L a I^ etc., which one would naively discard as a small part of the 
"free" Hamiltonian. That is, while total energy is conserved, the designation of an "interaction" part can 
be somewhat arbitrary, so it is best to be as consistent as possible and retain all terms in 1(4.3(1 . The 
electrostatic counterpart may be more familiar in the field, cf. the remark below 14.4(1 . 
For k t 1, 14.4(1 seems to diverge. However, physically this limit can only be achieved with two loops right 
on top of each other. The resulting /erromagnetic interaction is described by M < in ((4.2() . so that the 
last term in (|4.4|) tends to cancel the other two. Indeed, we have merely rewritten the finite l|4.3|l . 
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Two coupled 3JJ qubits can now be described by 

3 



22 



H 



EE 

c=a,b i=l 



Ql 

2C,;, 



— E ic cos ( 



M 



(4.5) 



with $ c = X]«=i 0ic/2e in ([4.4)1 . This description is valid for arbitrary inductances; the price 
to pay is that ()4.5|) involves six degrees of freedom. Let us turn to the quantum expansion 
right away, limiting ourselves to outlining the proper generalization of Section ILTTl 
In 0(L~ 1 ^ 2 ), the anisotropic oscillator problem is solved by 



-4> T A4> . 



(4.6) 



the inverse covariance matrix A is found from AC 1 A = L l /64e A , with C = diag(C a , C&) 
diagonal since the coupling between the qubits is purely inductive. One obtains 



A = -Lv'C [VChVC] - 1/2 VC , E zp = 4e 2 + ^) 
8e 2 V c a C b ) 



(4.7) 



Since even 2x2 matrix square roots are slightly tedious, we leave A unevaluated in most 
formulas [cf. below]. The correlation (entanglement) expressed by the Gaussian in ()4.6j) 
seems physically reasonable: for, e.g., again two small loops on top of each other, deviations 
of $ aj h from each other should be even more strongly suppressed than deviations from $ x . 
In £>(L°), one has 



. 2e 2 e ^ v x e - 2 ^A0 C -i v ^ e ^A^ 2 + e -<t>^ HQ _ Eo) ^ o = o 



(4.8) 



Since has two components, one needs a double integral jd(j) a d(j)be^^ TA ^ to cancel it, 
yielding a single solvability condition formally equal to ([3.8)1 . However, in the present case 
we additionally observe the decoupling 



Hq = H 0a + H( 



1 06 



Eq = E 0a + E { 



A)b 



= ^'oa(Xa, O a )ip' ob (xb, b ) . (4.9) 



Thus, the leading six-phase wave function tp can be strongly entangled, while the leading 
four-phase one ip' Q factorizes. Subsequently, ip 2 is obtained analogously to (|4.bjl . 
In (^(L 1 / 4 ), one again finds E\ — 0, while 



'03 



J 



.CA- 1 !^) 
16e 3 



(4.10) 



A Gaussian vector integral, followed by use of ()4.7j) . confirms the expected ((f)) = 2eLI^°\ 

In 0(L 1 / 2 ), the correction to H takes the form H 2a cf) 2 a + H 2h (f) 2 since (j4.5j) has no Josephson 
interaction [cf. the remark below (j4.6j) ]. This does not generate any interaction when <fi is 
integrated out with the entangled weight e~ 2< ^ TA< ^, so that 



i>2 = V4,(Xa, O a )lp 2b (Xb, Ob) , E 2 = E 2a + E 2b 

(H 0a -E 0a )^' 2a + [±(A- 1 ) aaJ £Z 2o -£ 2a ]V>oa = 



(4.11) 
(4.12) 



22 Not all of our assumptions (cf. note|2J) may hold in a figure-eight geometry |9( if part of the inductance is 
distributed along the shared leg. Therefore, this case deserves special attention. 
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and similarly for a «-> b. One could proceed to solve for ip4 but this will be omit- 
ted, since all that really matters is its relation to the effective wave function %p2 = 
/d0 a d0 fe ^4e-^//d0 a d0 b e- 2 ^. 

The next order 0(L 3 ^ 4 ) again is comparatively uninteresting, as it mainly justifies taking 



V* = o. 



Finally, in O(L), everything can be combined, and the effective Hamiltonian can be read 
off. Doing so leads to the central result of this paper, 



H = HL + HL 



'Oa ■ L1 0b 2^"* > ^ ' 

where H' 0a , H' ob are as in (|2.15jl (with XC = 2eA c B x ), except for the renormalizations 



Ei, 



E,„ = Ei, 



1 - 



Cc 

Cii 



2/A-l) 



(A 



a, ) 128 



1^2 

cc 



a, b) . 



(4.14) 



V. DISCUSSION 



One can view ()4.13|) in at least two ways. Experimentally, the bare couplings Ei C will 
usually be unknown and the diagonal contributions from the last term may be compara- 



tively small. In that case, most interesting will be the interaction MIa\xa, 9a) ^'(Xb, @b), 
derived with the antiferromagnetic sign without any remaining ambiguity. Also, this in- 
teraction is given in terms of explicitly known current operators, so that in e.g. the two- 
level (qubit) approximation on has H mt = M\I a \ \I b \ a^a b in the flux basis without addi- 
tional handwaving. Contrast 0, where improper averaging of the J^'s would have led to 

i.e., no interaction between the qubits. 
On the other hand, for detailed numerical comparison between the four- and six-phase 
theories, every last bit of parameter dependence needs to be made explicit. Therefore, we 
take the square root in ()4.7|) in the eigenbasis, yielding 



r(0), 



(A" 1 ), 



2y/2e 2 
C a 



I + LaC \ LbCb )v'L„C, l +L l ,C h +n 



+ 1 + 



n 

LbCb—L a C a 



LaCa+LfrCb—fl 



Q = y/(L a C a -L b C b ) 2 + 4M*C a C b , 
1 



det A 



C a Cb 

64? V L a L b (l—k 2 ) ' 



(5.1) 
(5.2) 

(5.3) 



while (A 

k 2 



_1 )b6 follows by a <->• b. The second square root in (|5.1|) has a positive argument 
since K" < 1 (positive-definiteness of L). Note that (|5.1|) reduces to A^ 1 as in (|3.1|l for 
Q — > ±(L a C a — L b Cb) (i.e., regardless of which of the uncoupled qubits has the larger LC). 

As already seen in Section 1111 CI the L-expansion involves two different dimensionless 
parameters — as expected in a system with Coulomb, Josephson, and magnetic energies. 
"Classical" self-flux effects go as LI 2 /Ei ~ (3. On the other hand, from (J3.26|) one reads off 
that the "quantum" zero-point effects go as e 2 ^ L/C ~ A -1 ~ Ec/ujlc- Thus, the former 
(latter) parameter does not contain the charging (Josephson) energy, which is also the reason 
why, to the given order, their contributions are separated in the final results. In particular, 
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for comparatively large junctions and not too small area, the O(L) self-flux effect could 
dominate over the 0(V~L) renormalization within the domain of validity of our expansion — 
as was implicitly used in Fig. ^ However, in the experiments for which all parameters 
are available, one has (LI c /<5> ,e 2 ^/L/C) = (0.003,0.015) H3 and (0.007,0.02) HQ 
respectively, so that, if anything, the renormalization effect dominates. 

Extension to > 2 qubits, using the same vector notation, is immediate. It is hoped 
that the methods presented will have additional use in other devices with near-constrained 
variables. 

I thank M.H.S. Amin (who also pointed out the exception in note I2*2*J) . M. Grajcar, J.F. 
Ralph, G. Rose, M.F.H. Steininger, M.C. Thorn, J.Q. You, and A.M. Zagoskin for fruitful 
discussions and comments on the manuscript. 
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